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ABSTRACT 


Similitude  theory  is  applied  to  the  problefn  of  fuel  sloshing  In 
accelerated  tanks  to  establish  criteria  for  the  design  of  model  eaperiments 
It  is  found  that  dynamic  modeling  is  possible  even  if  liquid  viscosity  is 
considered.  The  ranges  of  significant  parameters  and  the  selection  of 
model  liquids  are  discussed. 
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I.  INTRODUCTION 


Sloshing  motions  of  f'j«l  in  partislly  filled  tanks  may  established 
by  essentially  lateral  accelerations  cf  the  vehicle.  Such  fuel  motions  are 
of  importance  for  there  is  *he  possibility  of  extreme  osc.llations  if  the 
frequency  of  the  excitation  is  in  the  neighborhood  of  one  of  the  natural 
frequencies  of  the  liquid  fuel.  The  forces  exerted  by  the  fuel  on  the  tank 
walls  may  then  lead  tc-  pf  uroations  in  the  flight  trajectory,  or  may  even 
impose  severely  high  stresses  on  structural  components. 

Suppression  of  sloshing  modes  has  been  attempted  by  means  of 
various  mechanical  devices,  and, therefore, it  is  of  some  importance  to 
evaluate  the  effectiveness  of  ^uch  devices.  Purely  mathematical  approaches 
to  this  problem  are  exceedingly  difficult  because  of  the  boundary  conditions 
introdured  by  such  supprejsion  devices,  and  are,  therefore,  re  st  ricted 
cBsentidlly  to  very  simple  cases.  On  the  other  hand,  it  may  be  possible 
fa  evaluate  suppression  devices  by  me.in8  of  suitably  designed  and  con« 
ducted  model  tests. 

The  problem  of  concern  in  the  present  stuoy  is  the  specificalion  of 
criteria  for  model  tests  which  will  give  useful  quantitative  information 
ccacerning  sloshing  and  slosh  suppressoi  s  in  full-scale  fuel  tanks.  This 
is  to  be  accomplished  by  the  application  of  Similitude  Theory  to  obtain 


appropriate  expressions  for  the  model  -  prototype  relationships.  These 


may  in  turn  be  cmp'.oyed  to  yield  a  rational  deaign  of  model  teat*.  Thi*  ii 
in  contrast  to  previous  experimental  studies  of  the  fuel  sloshing  problem 
(Refs.  l»4),  which  have  apparently  been  conducted  without  regard  ftSr  simi* 
litude  requirements  and  consequently  are  believed  to.be  of  limited  useful* 
ness. 

The  purpose  of  the  present  report  is  to  distuti  the  basic  require* 
mertta  for  model  studies  of  fuel  sloshing  by  application  of  similitude  theory 


II.  SIMILITUDE  RELATIONS  FOR  FUEL  SLOSHING 


A.  General 

We  consider  the  vehicle  in  vertical  fltf:ht  with  constant  acceleration 
directed  alon^;  the  flight  path  (Fig.  I).  Sloshing  i*"  ’  stumed  to  be  esc ited 
by  tinne •dependent  accelerations  acting  essentially  normal  to  the  flight 
path.  For  purposes  of  the  present  study,  the  exciting  accelerations  will 
be  considered  either  as  purely  translational,  or  purely  rotational  about  .in 
axis  normal  to  the  flight  path. 

The  essential  features  of  the  sloshing  motion  are  further  governed 
by  the  physical  properties  of  the  wquid  and  the  pressurising  gas.  The 
importance  of  the  liquid  density  and  total  mass  are  obvious;  however, 
viscosity  and  surface  tensicn  require  .icditional  comment  If  suppression 
devices  are  incorporated  it  is  cle.ar  th.'t  their  effc  'tiveneas  miiut,  to  .* 
very  substantial  degree,  depend  upon  the  d.imping  f  .,rces  they  can  provide 
.IB  the  result  of  viscous  action  of  the  liquid.  If  models  of  very  small  sixe 
are  employed,  it  is  conct'ivable  that  the  effects  of  sloshing  may  be  aub« 
stantially  modified  by  surface  tension  forces  acting  on  the  suppression 
devices  and  at  the  tank  walls.  As  will  be  discussed  later,  it  does  not  seem 
possible  to  provide  rigorous  modeling  if  surface  tension  forces  are  included, 
and  further  it  is  believed  that  the  total  cont ribttion  of  such  forces  would  he 
small  compared  to  the  inertial  and  viscous  forces.  Consequently,  surface 
tension  forces  will  be  omitted  in  the  analysis. 
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Simiturlv.  II  iH  »uppoH«>d  th.it  th^  furcrt  produced  by  the  prcteurizin^ 
(•at  are  compared  ikith  the  other  forcec<  and  will, therefore, be 

omitted.  The  jj.ia  pressure  may  bo  accounted  for  implicUlv,  however.  «• 
it  is  inc hided  t.n  the  total  pressure  at  the  tank  wall  by  direct  enperpositton 
sinrr  the  jjas  volume  rein.imei  c  ssontiattv  constant  during  aloahing. 

Further  restrsrtians  to  the  analysts  art  provided  by  the  assumption 
of  smAll  exc I’ii'iyns.  rij;id  tank,  uiid  the  eicistence  of  geometrical  similarity 
in  all  ri  siiecla  between  ntotiel  ami  prototype.  The  suppression  devices 
ihf  niselves  will  he  et/ii^uiered  geametriCiiUy  sitiiilar  in  all  respects  between 
model  and  prototype,  ami  the  cof resjiondmg  inertial  and  apparent  mass 
forces  assuCiiilcU  with  them  will  be  neglected  in  comparison  with  the  vis-  . 
coiis  forces. 

B.  Simtlttudo Jleiatuins  -  Trartslatie>.-.al  Kxeit.ition 

TUv  relevar.t  parameters,  basei^on^the  p/.^eedjog  discussion,  .will.  . 
be  taken  to  bo  (Fig.  1).  , 

a  longitudifi.ll  acceleration  acting  on  tank  (LT"^ 
d  ~  tank  djatneter  f  1*1 

f  s  resultant  liquid  force  on  tank  wall  (MLT”^) 
h  ^  depth  of  liquid  in  tank  (L) 

Xj,  =•  excitation  amplitude  (Li 
p  -  liquid  viscoady  fkil,"  *T‘*) 

[>  =  liquid  density 

r  5^  excitation  period  iTj 


i 

I 


five  dimcnsionlcfls  groups.  For  c.it  h  r*:'  thcH<*  uroup*  fn-r  of  fh**  tirKnowi-na 
ntust  he  set  arbitrarily.  Since  the  rcKultani  liquid  force  on  the  tank  wall 
is  the  d.'pcndcnt  ••anable  of  the  problem,  we  spiy  assign  >  the  value  of 
unity  in  one  solution  and  zero  in  all  olh<:r8,  thus  insuring  that  F  occurs 


only  once  and  only  to  the  firs*  power. 
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Number. 


Should  interest  reside  not  in  the  rosuiiant  lii^utd  force  on  the  i.mk 
wall,  but  rather  in  the  resultant  liquid  prei«»ure.  a  similar  anslyai*  «ould 
yield 


so  that 

n.  ,  Ell  =  — li — 

*  p<i^  p(d/T)^ 


(12) 


which  is  in  the  form  of  a  pressure  coefficient  (Euler  Number). 

The  preceeding  analysis  has  considered  the  excitation  to  ’te  a  dls> 
placement,  in  anticipation  of  actual  model  experiments.  The  exctiation 
could,  however,  have  been  taken  equally  well  as  a  force  ^  (or  even  an 
acceleration),  in  which  case  the  appru,)tiatc  dime  sionless  proup  would 
become 


P 

(pdh  (d/T^) 


(13) 


Consideratiun  of  the  surface  tension  in  the  foregoing  analysis 
would  have  resulted  in  a  new  IT  group  of  the  formpd'/er^,  which  corre¬ 
sponds  to  a  Weber  Number.  As  will  be  seen  from  later  considerations, 
the  simultaneous  solution  of  this  new  group  and  11^  would  be  exceedingly 
difficult,  in  view  of  the  model  liquids  readily  available.  Further,  we 
assume  that  surface  tension  forces  are  sm.ill  iomp.ired  with  the  inertisl 
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and  viHcouii  force*,  and  conarquently  wr  *hall  omit  any  additional  considera¬ 
tion  of  ■iirfaco  tension. 

C .  Similitude  Relaticnw  -  Rotational  Eacitation 

In  till*  cA*r  the  pertinent  parameters  remain  the  same  as, in  the 
tran-tlatinnal  case  except  that  the  excitation  is  how  defined  by  two  para¬ 
meter*.  the  anKul.ir  rotaticii  and  the  location  of  the  totational  axis  b 
{Fi^.  Jl.  Since  u-e  consider  only  small  excitation  amplitudes,  it  may  be 
noted  that  ^  and.  iherefure.the  general  functional  Equation  (10) 
applies  to  both  kinds  of  excitation.  The  use  of  a  moment  M  *  Pb  a*  the 
excitation  would  result  in  a  group  formed  by  fl^  multiplied  by  the  ratio 


b/d. 
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III.  MODELING  CONSIDERATIONS  BASED  ON  THE 
SIMILITUDE  RELATIONS 

A.  Grnv*ral 

Equation  (10),  and  the  associated  alternate  expressions  for  special 
canes,  governs  the  design  of  models  to  simulate  a  given  prototype.  If  all 
dimensionless  groups  in  this  equation  have  the  same  numerical  value  for 
both  model  and  prototype,  then  the  forces  and  pressures  measured  on  the 
model  are  directly  applicable  to  the  prototype. 

Denoting  the  ratio  of  model  to  prototype  parameters  by  the  sub* 
script  r,  we  note  that  fluid  depth  and  excitation  displacement  amplitude 
are  scaled  in  the  name  proportion  ar  the  geometrical  scale.  Thus 
dr*dm/‘>p 

and 

I 

hr  =  ^‘or  •  br  *  df  ,  0,  *  I  (15) 

The  remaining  model  parameters  are  to  be  established  by  further 
study  ol  Eq.  ( 10) . 

B .  Modeling  Considerations  -  Viscosity  Neglected 

It  IS  instructive  to  consider  the  modeling  parameters  in  the  case 
fh,at  VISCOUS  effects  are  neglected  throughout.  This  means  that  Ilj  drops 

out  of  the  previous  analysis,  and  the  relation  between  model  and  prototype 
depends  on  the  geometric  factors  discussed  above  and  the  time  acale  factor 

tJ  *  d,/a. 


(16) 
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This  result  means  that  any  sire  model  in  a  accrler.ttton  field  ean  he 
used  to  simulate  any  sise  prototype  in  anv  acceleration  field  by  varying 
the  time  scale,  and  would  further  permit  the  use  of  any  model  liquid  since 
the  density  appears  only  in  the  force  (or  press:ire)  parameter. 

In  presenting  measured  data,  the  parameter  ar^/d  may  be  used 
as  the  primary  independent  variable.  Since  frequency  is  a  more  com* 
monly  used  variable,  it  may  be  introduced  lor  the  period  to  give  dw^/a. 
The  force  (or  pressure)  group  may  also  be  modified  by  introduction  of 


£q,  ( 16)  to  give 


apd 


T 


(in 


the  measured  data  may  then  be  presented  with  F/apd^  and  dw^/a  as  the 
primary  variables  and  the  remaining  quantities  as  parameters, 

C .  Modeling  Considerations  -  Viscosity  Included 

With  viscous  forces  ircluded,  the  data  nie>  still  be  presented  aa 
indicated  above;  however,  the  model  design  is  now  subject  to  certain 
restrictions  arising  from  n|^,  the  Reynold's  Number.  Introducing  Eq.  (lb). 


we  find 


from  which 


•^r  • 


'A 


-1/3 


f  IB) 


(19) 


2/3 
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ThiM  imporlant  rrttuU  iihoH-!i  that  the  geometrical  scale  is  determined  by 
the  acceleration,  liquid  density,  and  liquid  viscosity  ratios.  If  it  is 

.d 

desired  to  test  models  so  as  to  cover  a  range  of  prototype  accelerations, 
we  m.iy  either  fix  the  model  diameter  and  vary  the  liquid  properties,  or 
fix  the  model  liquid  and  vary  the  model  sise.  Since  the  range  of  liquid 
properties  is  very  restrict!  d.  the  latter  course  is  chosen. 

s 

In  order  to  employ  models  of  reasonably  small  sise  to  study  large 
prototviM's  ins|ii*clion  of  Er  ( |0)  indic.ates  that  the  model  liquids  must  have 
high  density  or  low  viscosity,  or  both,  compared  to  the  prototype  liquid. 
Mercury,  of  course,  has  a  very  high  density,  but  is  undesirable  for  a 
variety  of  reasons;  some  of  the  organic  solvents  have  very  low  viscosity, 
therefore.offer  tome  promise. 

As  an  example  of  model  simulation,  assume  the  following  typical 

case: 


Parameter  Prototype 

Liquid  Kerosene 

Specific  gravity  .83 

Viscosity  2.SC.P. 

Acceleration  l-5g 

Diameter  9  ft 


For  a  prototype  acceler.ation  of  1g, 


1 

4.21 


Model 

Methylene  Chloride 
1.336 
.465  c.p. 


I ) 


ThcTtforc,  a  9-ft  diamrlor  prototyp«*  with  kcropcnc  futl  and  ttndrritoing 
acrclrrationk  from  l-Sg  may  bo  simulated  by  a  aoru  a  o*  modola  of  Z  14 
to  3.66-ft  diam'^tcr  with  mothylcno  chloride  as  th'  model  liquid.  The 
model  frequonry  rAr.;;o  la  determined  from  E^.  116)  as 
=  1/t^  =  (ar/dp)*^^  *  2.052  to  0.701 

It  ia  clear  from  the  foregoing  that  the  inclusion  of  an  additional 
group  to  account  for  surface  tension  would  have  rendered  the  model  design 
virtually  impossible  because  of  the  inability  to  find  a  model  liquid  of  appro* 
priate  properties  to  satisfy  «.Il  dimensio.-iless  (.roups  simultaneoualy,  in  a 
small  model. 

It  may  be  noted  in  passing  that  the  ekperimenlt  reported  m  Refs. 

1*3  were  restricted  to  small  models  with  w’ater  as  the  model  liquid  n.isrd 
on  the  considerations  of  t!i«.  present  report,  t.hose  reuulla  cannot  adequately 
simulate  sloshing  with  supjiressors  because  of  the  f;«ilurc  to  model  mscous 
forces.  While  liquids  other  than  water  were  also  employed  in  lh«  expert* 
ments  of  Ref.  4,  the  liquid  properties  were  not  properly  fitted  to  the 
model  sires  and  hence  the  results  cannot  be  usid  to  ev.iluate  supjiression 


devices. 
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